Truncated-Bloch-wave solitons in optical lattices 
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We study self-trapped localized nonlinear states in the form of truncated Bloch waves in one- 
dimensional optical lattices, which appear in the gaps of the linear bandgap spectrum. We demon- 
strate the existence of families of such localized states which differ by the number of intensity peaks. 
These families do not bifurcate from the band edge, and their power curves exhibit double branches. 
Linear stability analysis demonstrates that in deep lattice potentials the states corresponding to 
the lower branches are stable, whereas those corresponding to the upper branches are unstable, 
independently of the number of peaks. 
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The theoretical study of localized states in periodic po- 
tentials is driven by many experimental demonstrations 
in nonlinear optics [ll and the physics of Bose-Einstein 
condensates (BECs) [2, i, S H- One of the most strik- 
ing effects is the possibility of spatial localization due to 
an interplay between periodicity and nonlinearity, even 
when the nonlinearity is repulsive. The corresponding lo- 
calized states have long been known in nonlinear optics as 
gap solitons [6| whose frequencies lie inside the gaps of the 
linear band-gap spectrum of a periodic structure. Simi- 
lar objects have been observed in atomic physics, where 
the experimental generation of BEC gap solitons proved 
to be a challenge due to the requirement of low atom 
numbers and densities The physics of such localized 
states is well understood, and usually the gap solitons 
appear through bifurcations of the localized states from 
the edges of the bands of periodic solutions Q ■ 

Seemingly different localized states with steep edges 
and a large number of atoms populated evenly in a num- 
ber of adjacent lattice sites was observed a few years ago 
in BEC loaded into an optical lattice [|[ . It was observed 
that in a deep one-dimensional optical lattice the BEC 
wave packet does not diffuse, but instead its initial ex- 
pansion stops and the width remains finite. This effect 
has been attributed to the self-trapping mechanism of en- 
ergy localization and the existence of novel types of broad 
nonlinear localized states, in the form of truncated non- 
linear Bloch waves in the gaps of the matter- wave linear 
spectrum 0] . In particular, Alexander et ai. 0] revealed 
that these robust nonlinear localized states may exist in 
all dimensions and can have arbitrary extension within 
the lattice. They seem to provide an important missing 
link between gap solitons and nonlinear Bloch waves, and 
can be termed truncated- Bloch-wave solitons. 

In this paper, we consider a model for the Bose- 
Einstein condensate in a one-dimensional optical lattice 
under repulsive nonlinearity, and systematically study 
the existence and stability of the truncated-Bloch-wave 
solitons. On the question of existence, we reveal that 
these solitons exist as countable families which are char- 
acterized by the number of intensity peaks. These so- 



lution families do not bifurcate from the edges of Bloch 
bands, and their power curves exhibit double branches. 
If the depth of the lattice potential falls below a cer- 
tain threshold value, these truncated-Bloch-wave solitons 
cease to exist. We also demonstrate that in a deep lat- 
tice potential these solitons on the lower branches of 
the power curves are stable, while those on the upper 
branches of the power curves are unstable, independent 
of the number of intensity peaks. In a shallow lattice 
potential, both branches of these solitons are linearly 
unstable. These results provide a solid background for 
the earlier experimental observations of broad localized 
states in BECs Q as well as the corresponding numerical 
studies @, and they reveal the specific role these novel 
types of self-trapped localized states play as a link be- 
tween nonlinear Bloch waves and gap solitons in periodic 
potentials. 

The physical system we consider is that of a Bose- 
Einstein condensate loaded into a one-dimensional opti- 
cal lattice under repulsive nonlinearity 0, [1[ . The math- 
ematical model for this system is the one-dimensional 
Gross-Pitaevskii equation with a lattice potential [l(| 



iU t + U xx - (V sin 2 x)U - \U\ 2 U = 0, 



(1) 



where U is the mean- field wavefunction of the BEC, and 
the lattice potential is 7r-periodic with potential depth Vo- 
We have used the same normalization as that found in 
Ref. 0] to make the model dimensionless, with all quanti- 
ties in units of characteristic scales of the lattice. In par- 



ticular x and t are in units of = d/ir and < 



h/E r 



respectively, where d is the period of the lattice and 
E rec = h 2 /2ma 2 L is the recoil energy for an atom of mass 
m absorbing a lattice photon. The depth of the lattice, 
Vo, is in units of E rec . We have assumed that the con- 
densate is tightly confined in the transverse dimensions 
to use the standard dimensionality reduction procedure 
to reduce the dynamics to an effectively one-dimensional 
model. 

Solitons in this system are sought in the form of 
U(x,t) = u(x)e~ lfJ,t , where fi is the chemical potential, 
and the spatial function u(x) is real-valued and satisfies 
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FIG. 1: Nonlinear Bloch waves bifurcated from the lower edge 
of the first Bloch band under repulsive nonlinearity (Vo = 6). 
(a) The amplitude versus the chemical potential diagram; (b, 
c) nonlinear Bloch waves at two chemical potentials marked 
in (a). The shaded stripes represent low-potential regions. 

the reduced equation 

Uxx — (Vo sin 2 x)u — u 3 = —/iu. (2) 

In this paper, we first consider a deep lattice potential 
with Vo = 6, where the existence and linear stability 
of truncatcd-Bloch-wave solitons will be determined in 
detail. These solitons and their stability in shallower po- 
tentials will be addressed toward the end of the paper. 

To understand the origin of the truncated-Bloch-wave 
solitons, we first examine periodic nonlinear Bloch waves 
which originate from linear Bloch waves at the band 
edges. In a deep potential with Vo = 6, the linear dis- 
persion curves of Eq. @ can be found in Fig. HJa) of 
Ref. [ll| . At the lower edge of the first Bloch band [which 
is 2.0632 in Fig. UK a)], Eq. ([2]) admits an infinitesimal 
(linear) Bloch wave which is 7r-periodic, and its adjacent 
intensity peaks are in-phase with each other. This lin- 
ear Bloch wave (in arbitrary units) is shown in Fig. 2(1) 
of Ref. [ll| and appears similar to the profile shown in 
Fig. [TJb) here. When increases from this band edge 
(into the first Bloch band), this linear Bloch wave bi- 
furcates into a nonlinear 7r-pcriodic Bloch wave with a 
finite amplitude (under the present repulsive nonlinear- 
ity). Numerically we can easily compute these nonlinear 
Bloch waves, and their amplitudes versus the chemical 
potential /i are displayed in Fig. [lja). Here the ampli- 
tude of a nonlinear Bloch wave u(x) is defined as the 
maximum value of u(x). At two representative chemical 
potentials fi = 2.16 and 4 (the former lies inside the first 
Bloch band, and the latter lies inside the first bandgap), 
the corresponding nonlinear Bloch waves are displayed in 
Figs. [ljb,c), respectively. As can be seen, as fi increases 
from the lower band edge, the amplitude of the nonlinear 
Bloch wave also increases. 

Now we truncate these infinitely extended nonlinear 
Bloch waves to a finite number of intensity peaks (with 
the [i value fixed), and ask whether the truncated Bloch 
waves can develop exponentially decaying tails at large 
\x\ values so that they become localized soliton solutions 
of Eq. ||2J). When the chemical potential lies inside the 
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FIG. 2: Truncated-Bloch-wave solitons on the upper and 
lower power branches of the three-peak (left column) and 
seven-peak (right column) families in the first bandgap (Vo = 
6). The chemical potentials of these solitons are marked 
in Fig. [3] In (d), the corresponding nonlinear Bloch wave 
(dashed curve) is also shown for comparison [see also Fig. 

arc)]. 

first Bloch band [see Fig. [31b)], the truncated Bloch 
waves will still possess continuous- wave tails and are thus 
still infinitely extended. But when the chemical poten- 
tial lies inside the first band gap [see Fig. Etc)], then due 
to repeated Bragg reflections, the truncated Bloch waves 
can decay exponentially at large |x| values and thus be- 
come localized solitons. To confirm this, we numerically 
compute soliton solutions of Eq. ([2]) by the sq uared- 
operator iteration methods developed in Ref. [l2j using 
the truncated nonlinear Bloch waves as initial conditions. 
When truncating the Bloch wave of Fig. [ljc) to three and 
seven peaks, we indeed obtained soliton solutions which 
are displayed in Fig. [Hb,d) respectively. These solitons 
closely resemble the corresponding nonlinear Bloch wave 
[see Fig. rjjd) for comparison], but their tails arc now 
exponentially decaying rather than infinitely extended. 
If the Bloch wave of Fig. [IJc) is truncated to any other 
number of peaks, the corresponding soliton state could 
be obtained as well. Thus a countable number of solitons 
can be found from different truncations of a nonlinear 
Bloch wave (at the same \i value). These solitons are 
the ones which we termed truncated-Bloch-wave solitons 
above. They were first observed in BEC experiments Q 
and then found numerically in Ref. Q. 

When the chemical potential changes, the truncated- 
Bloch-wave solitons will deform. Hence each truncated- 
Bloch-wave soliton generates a continuous family of such 
solitons parameterized by the chemical potential. We 
will call the solution family of the three-peak truncated- 
Bloch-wave soliton in Fig. (Hb) the three-peak family, 
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and that of the seven-peak soliton in Fig. [2jd) the seven- 
peak family, and so on. Important questions which have 
not been addressed before arc: what solutions are con- 
tained in each family? and what do their power dia- 
grams look like? Here the power of a soliton is defined 
as P = \u\ 2 dx as usual. To address these ques- 
tions, we have numerically obtained the entire families 
of two-peak to seven-peak solitons. The power diagrams 
of three-peak and seven-peak families are displayed in 
Fig. El (power curves for other-peak families are simi- 
lar). These power diagrams show three important com- 
mon features about truncated-Bloch-wave solitons. The 
first feature is that these solitons exist only inside the 
bandgap, but not inside a Bloch band. This contrasts 
with the associated nonlinear Bloch waves which exist in 
both the bandgap and the Bloch band [see Fig. da)]. 
The second feature is that these solitons do not bifurcate 
from the (upper) edge of the (first) Bloch band, thus 
they are fundamentally different from the gap solitons 
studied before Q which do bifurcate from edges of Bloch 
bands. The reason these solitons can not bifurcate from 
the upper edge of the first Bloch band is obviously that 
adjacent peaks in these solitons are in-phase, which does 
not match the Bloch wave at the band edge (where ad- 
jacent peaks arc out-of-phasc). Note that these solitons 
are closely related to nonlinear Bloch waves which bi- 
furcate from the lower edge of the first Bloch band (see 
Fig. []}. But unlike those nonlinear Bloch waves, these 
solitons can not exist inside the first Bloch band, thus 
they can not bifurcate from the lower edge of the first 
Bloch band. The third common feature of the truncated 
Bloch-wave solitons is that their power diagrams exhibit 
double branches. This double-branch phenomenon has 
been reported for dipole and vortex solitons in periodic 
potentials recently [HI, [l4| ■ Now we know that this phe- 
nomenon will always occur for any soliton family which 
does not bifurcate from edges of Bloch bands (or edges 
of the continuous spectrum) as is the case for truncated- 
Bloch-wave solitons. The three-peak and seven-peak soli- 
tons shown in Fig. [2jb,d) lie on the lower branches of 
their solution families (more specifically at points 'b, d' 
marked in Fig. [3]). On the upper branches, the solitons 
will develop two extra intensity peaks at the two edges of 
the solitons, and these extra peaks are out-of-phase with 
the other intensity peaks in the interior. Two examples 
on the upper branches of the three-peak and seven-peak 
families are displayed in Fig. IHa,c), and their positions 
on the power diagrams are marked in Fig. [3J It is noted 
that these solution families of various peaks are distinct 
from each other, with no s-shaped bifurcations connect- 
ing them near the band edge in the power diagram. When 
the chemical potential fi approaches the right end of the 
first bandgap, all solutions become less localized. If \i 
crosses this right end into the (second) Bloch band, solu- 
tions become delocalizcd due to resonance with the Bloch 
modes, and the powers of the solutions become infinite. 

Stability of these truncatcd-Bloch-wave solitons is 
clearly an important issue, and this issue has not been 
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FIG. 3: Power curves of the three-peak (a) and seven-peak 
(b) families of truncated-Bloch-wave solitons (Vb = 6). The 
solutions at marked points are displayed in Fig. [2] 



considered before in the literature. Here we investigate 
the linear stability of these solitons. Following the 
standard procedure, we perturb the soliton as U(x,t) = 
e-w* [u(x) + (v(x) - w{x))e xt + (v*(x) + w*(x))e x ' t ] , 
where v, w <C 1 arc normal-mode perturbations, A is 
the eigenvalue of the normal mode, and the superscript 
'*' represents complex conjugation. Inserting this 
perturbed soliton into the evolution equation {T]) and 
dropping higher-order terms in (v, w), a linear eigenvalue 
problem L(v,w) T = X(v,w) T is obtained. Here L is the 
linearization operator, and the superscript 'T' represents 
the transpose of a vector. Linear stability of a soliton is 
determined by the spectrum of the linearization operator 
L, and the existence of any eigenvalue with a positive 
real part implies linear instability of the soliton. Wc 
have computed this spectrum for truncated-Bloch-wave 
solitons of various families in the deep potential of 
Vo = 6 (see Figs. [5] and [3J by the Fourier collocation 
method [15j , We have found that for all these solution 
families, solitons on the lower branches are always stable, 
and those on the upper branches always unstable. To 
illustrate, we pick the two solitons in Fig.[2)Jc, d), which 
lie on the upper and lower branches of the seven-peak 
family in Fig. [3Jb) . The linear-stability spectra of these 
two solitons are shown in Fig. BJb,a), respectively. We 
see that eigenvalues of the lower-branch soliton all lie 
on the imaginary axis, indicating that this soliton is 
linearly stable. But the spectrum for the upper-branch 
soliton contains real positive eigenvalues, indicating 
its linear instability. Intuitively, the linear instability 
of this upper-branch soliton is easy to understand. 
Indeed, this soliton has two intensity peaks on the two 
edges which arc out-of-phase with the interior peaks. 
Thus the structures near the edges of this soliton can 
be viewed as out-of-phase dipoles. It is well known 
that in lattices under attractive nonlinearity, in-phase 
dipoles are linearly unstable (see, e.g., Ref. [Hi). If the 
nonlinearity is repulsive as it is in the model ([1]), dipole 
stability is switched, and out-of-phase dipoles become 
linearly unstable (see e.g. Ref. [17|). Consequently, this 
upper-branch soliton in Fig. [He) should be linearly un- 
stable. Regarding the linear stability of the lower-branch 
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FIG. 4: Linear-stability analysis of the seven-peak family of 
truncated-Bloch-wave solitons in Figs. 1 and 2 (Vb = 6): (a) 
stability spectrum of the lower-branch soliton shown in Fig. 
[2jd); (b) stability spectrum of the upper-branch soliton shown 
in Fig. Etc). 



soliton, however, it is less obvious intuitively, and this 
result is one of the main findings of this paper. 

The above linear-stability result of lower-branch 
truncated-Bloch-wave solitons in a deep lattice poten- 
tial shows that these solitons should also be nonlinearly 
stable under weak perturbations. To confirm this, we 
take the seven-peak soliton of Fig. HJd), and modulate 
its seven peak amplitudes by 10% random perturba- 
tions. Specifically we add to this soliton a perturbation 
which is a superposition of seven Gaussian humps of the 
shape exp(— x 2 ) centered at the seven lattice sites of the 
soliton. The amplitudes of these Gaussian humps are 
taken randomly at the level of about 10% of the soli- 
ton's amplitude. For the realization of amplitude values 
(0,-0.08,0.06,-0.06,0.01,0.05,-0.11) for these Gaus- 
sian humps, the nonlinear evolution of the perturbed soli- 
ton after 60 time-unit evolution is displayed in Fig. [5] We 
see that the evolution is stable as expected, showing the 
robustness of these new types of soliton structures in lat- 
tice potentials. Evolution with other realizations of the 
perturbations is qualitatively similar. 

These theoretical results appear very similar to the 
experimental observation of nonlinear self-trapping of 
Bose-condensed 87 Rb atoms in a one-dimensional deep 
periodic potential @. In the experiment of Ref. 
it was observed that with a higher number of atoms, 
the repulsive atom-atom interaction led to the forma- 
tion of a self-trapped broad localized state with steep 
edges and roughly uniform interior, closely resembling 
the truncated-Bloch-wave soliton shown in Fig. [UJd) . We 
note however that the lattice depth used in the exper- 
iment corresponds to Vq ~ 11, much deeper than the 
Vb = 6 considered in Fig. O Furthermore, as noted 
in Ref. [|[, the experiment was not carried out in the 
one-dimensional regime, with higher-order transverse ef- 
fects playing an important role. To directly compare 
our results with those of experiment we may obtain the 
physical density and atom number from our normal- 
ized model by multiplying the normalized density \U\ 2 
by 2 x 10~ 6 /87rafja s |, where a s is the s-wave scattering 
length, to obtain the density in atoms. cm -3 and multi- 
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FIG. 5: Nonlinear evolution of the seven-peak truncated- 
Bloch-wave soliton of Fig. [2jd) under 10% perturbations: 
(a) intensity (|w| 2 ) contours in the (x,t) plane; (b) initial per- 
turbed state of |u|; (c) final state of \u\ after 60 time-unit 
evolution. 



ply the power P by H/(4:\a s \muj±aL), where uj± is the 
transverse trap frequency, to obtain the total number 
of atoms. Inherent in our one-dimensional approxima- 
tion is the assumption that the nonlinear energy at the 
peak density is much less than the transverse trap en- 
ergy, i.e. that |u pea fc| 2 ft/(TOa| / u;j^) << 1, and as such 
we consider the experimental parameters d = 3.5/im and 
lu± = 320Hz. For a 87 Rb BEC the seven-peak solution 
shown in Fig. Hfd) then contains approximately 400 
atoms, with ~ 60 atoms per lattice period and a peak 
density of ~ 2.8 x 10 13 atoms.cm~ 3 , all experimentally 
reasonable values. 

The above theoretical results pertain to truncated- 
Bloch-wave solitons in a deep lattice potential. What will 
happen to these solitons if the lattice potential becomes 
shallow? To address this question, we take a shallow po- 
tential with Vo = 1, and investigate the existence and 
stability properties of the truncated-Bloch-wave solitons. 
For the seven-peak family, the results are summarized in 
Fig. [6] (results for other-peak families are qualitatively 
similar). In this shallow potential, we see from Fig.^a) 
that these solitons exist in a smaller region of the first 
gap. From Fig. [J5{b,c) we see that the edges of these 
solitons become less steep, their tails become longer, and 
solution profiles on the upper and lower branches be- 
come more similar. From Fig. E^d), we see that the 
soliton on the lower branch now becomes linearly un- 
stable due to oscillatory instabilities caused by complex 
unstable eigenvalues (this holds for other solitons on the 
lower branch as well). The solitons on the upper branch 
are certainly also unstable. Thus in a shallow potential, 
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FIG. 6: Truncated-Bloch-wave solitons in a shallow lattice 
potential (Vb = 1): (a) the power diagram; (b, c) solitons on 
the upper and lower branches as marked in (a); (d) linear- 
stability spectrum for the lower- branch soliton in (c). 

both branches of truncated-Bloch-wave solitons are lin- 
early unstable, which contrasts with the deep-potential 
case. 

By comparing truncated-Bloch-wave solitons in the 
deep and shallow potentials of Vb = 6 and Vq = 1, we see 
that when the potential becomes shallow, the existence 
region of these solitons shrinks toward the upper part of 
the (first) bandgap (see Figs. [3] and 0. This shrinkage 
continues until the potential depth drops below a criti- 
cal value of Vb = 0.65, where these truncated-Bloch-wave 
solitons all disappear (see also Ref. @). Thus very shal- 
low lattice potentials can not support truncated-Bloch- 
wave solitons. Regarding the linear stability of these soli- 
tons, the results above show that in deep lattice poten- 
tials, lower branches of truncated-Bloch-wave solitons are 
stable, but in shallow potentials, both branches of these 
solitons are unstable. As the potential depth decreases, 
the stable region of truncated-Bloch-wave solitons grad- 
ually shrinks toward the middle part of the lower branch, 
and then vanishes when Vq ~ 2.3. Thus whole families of 
these solitons become unstable when Vq drops below ap- 
proximately 2.3. The experiment of Ref. @, which used 
Vb ~ 11, was thus far above this instability threshold. 



As Vb is defined in the unit of the lattice recoil energy 
E rec , the dimensionless instability threshold of Vb ~ 2.3 
thus corresponds to the dimensional instability threshold 
of Vb ~ ll.35h 2 /md 2 , which is inversely proportional to 
the atom mass and the square of the lattice period. 

In this paper, we only analyzed truncated-Bloch-wave 
solitons in the first bandgap under repulsive nonlincar- 
ity These solitons were obtained by truncating the non- 
linear Bloch waves originating from the lower band edge 
of the first Bloch band. It should be noted that this 
branch of nonlinear Bloch waves (see Fig. [I} contin- 
ues into higher Bloch bands and bandgaps. In higher 
bandgaps, by truncating these nonlinear Bloch waves, we 
may obtain higher-gap truncated-Bloch-wave solitons. It 
is also noted that from every edge of a Bloch band, a 
branch of nonlinear Bloch waves bifurcates out for both 
attractive and repulsive nonlincaritics [the bifurcation 
goes rightward (see Fig. [Ha)) under repulsive nonlincar- 
ity and leftward under attractive nonlinearity] . When- 
ever a branch of nonlinear Bloch waves enters a bandgap, 
we may obtain the corresponding truncated-Bloch-wave 
solitons. These other truncated-Bloch-wave solitons can 
be similarly analyzed, but this lies outside the scope of 
the present paper. 

In conclusion, we have analyzed the existence and sta- 
bility of truncated-Bloch-wave solitons existing in the 
gaps of the linear band-gap spectra of matter waves. We 
have demonstrated that such self-trapped structures form 
families of broad but localized solutions which can be 
identified by the number of peaks. These localized states 
do not bifurcate from the band edge, and their power di- 
agrams exhibit double branches. Linear stability analysis 
demonstrates that in a deep lattice potential, the solu- 
tions on the lower branches are stable regardless of the 
number of peaks; but in a shallow potential, solutions 
on both branches are unstable. When the lattice poten- 
tial becomes weak enough, these truncated-Bloch-wave 
solitons cease to exist. 
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